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ABSTRACT : 

The production of ultrahigh energy photons, electrons 
and neutrinos as the decay products of pions produced in 
photomeson interactions between cosmic-ray nucleons and the 
blackbody microwave background is discussed in terms of the 
resultant energy spectra of these particles. Simple 
asymptotic formulas are given for calculating the ultrahigh 
energy photon spectrum predicted for the universal cosmic-ray 
hypothesis and the resulting spectra are compared with 
those obtained previously by numerical means using a 
different propagation equation for the photons. Approximate 
analytic solutions for the photon spectra are given in terms 
of simple power-law energy functions and slowly varying 
logarithmic functions. These can be used to estimate ultra- 
high energy photon fluxes for various astrophysical parameters. 
The generic relation between the various secondary components 
is then discussed in terms of their astrophysical implications 
which are summarized in the conclusion. 
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The question of the origin of ultrahigh energy cosmic- 
rays evidenced by massive air-showers is still a matter of 
controversy. The older view of an extragalactic origin has 
recently been summarized by Brecher and Burbidge (1972) . 

A more recent view of galactic origin has been argued for by 
Syrovatsky (1971) and by Stecker (1971) , and we refer the 
reader to these references for the pertinent background 
discussion. The purpose of this paper is to discuss a 
possible test of the universal origin for ultrahigh energy 
cosmic-rays and also for the neutrino origin hypothesis 
of Berezinsky and Zatsepin (1969 a,b ) . This test lies 

in a search for the ultrahigh energy photons which must be 
produced by photomeson interactions between "universal" 
ultrahigh energy cosmic-rays and the universal 2.7K microwave 
background radiation if both are present throughout the 
universe. 

It was pointed out by Gould and Schreder (1966, 1967) 
and Jelly (1966) that photons with energy above 10 14 eV 
will be seriously attenuated by pair-production interactions 
with the 2.7K background. However, recently Bonometto 


(1971) has pointed out that at the universe may be some two 

/ 
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orders of magnitude more transparent to 10 1 photons than 
had previously been estimated due to the development of an 
electromagnetic cascade process in the 2. IK radiation medium. 

In this cascade process, the initial photon interacts with a 
2.7K background photon to produce an electron-positron pair, 
giving most of its energy to one of the electrons. This 
electron, in turn, Compton scatters with a 2.7K background 
photon to produce another photon with almost the same energy 
as the electron. Thus, the total effect of the two interactions 
is the same as if the original photon had been somewhat 
degraded in energy rather than completely absorbed. Thus, 
as Bonometto points out, this allows the transmission of 
photons of energy E^,^10 19 eV from effective distances of over 
^ 4xl0 26 cm (120Mpc). (Allcock and Wdowczyk (1972) considered 
this process in more detail). As a consequence, Wdowczyk, et al . 
(1971 a,b) have pointed out the possibility of observing 
photon initiated showers in good proportion to the total number 
of showers at energies ^10 19 eV. It thus becomes important 
to estimate the intensity of these photons as Wdowczyk, et al. 
have done and to consider their generic relation to other 
possible components of ultrahigh energy cosmic radiation such 
as the neutrinos considered by Berezinsky and Zatsepin (1969a, b) 
whose mean-free-path is larger than the visible radius of the 
universe. In this paper, estimates are made of the 
ultrahigh energy photon, electron, and neutrino fluxes with 
assumptions which differ somewhat from those of Wdowczyk, 
et al. Simple asymptotic formulas are given for calculating 
the ultrahigh energy photon spectrum taking cascading into 
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account and the results are compared with those of Wdowczyk, 
et al. (1971 a,b; 1972) and discussed. The approximate 
analytic solutions given here in terms of simple power-law 
energy functions and slowly varying logarithmic, functions 
can be used to estimate fluxes of ultrahigh energy photons 
and other secondaries for various astrophysical parameters 
without resorting to numerical solutions. 

The main source of ultrahigh energy photons is photo- 
meson production interactions between ultrahigh energy cosmic- 
rays and blackbody photons of the 2.7K microwave background 
radiation i.e. interactions of the type y + p->TT°+p. This 
source may produce a detectable flux of 10 19 - 10 20 eV 
photons if both the 2.7K background and the ultrahigh 
energy cosmic rays are universal. To derive an approximate 
source function for the production of y-rays in photomeson 
interactions with microwave blackbody photons, we assume 
all the photons to be at the average energy £ q ~ 2.7kT - 
6.4x10 4 eV so that 

( 1 ) 

W £) = n bb 6(e_e o ) ' n bb " 400cm 3 

The energy of the photon in the cosmic-ray proton rest 
system is 

(E p /Mp)e o (l-cos0 ) 


e 


( 2 ) 



4 


There is a large peak in the photomeson production cross 
section at £* = (). 35M p due to the A (1238) resonance and 
since the cosmic-ray spectrum drops off rapidly with increasing 
energy, most of the pion production occurs at this resonance 
energy. Thus, from equation (2) we make the approximation 


a (c') - 


■ 

x- 


0.35M 
l 

e E 
o p 


o o s<x -VV 


(3) 


where 


x = 1-COS0, E E Q : 35 ^ -- -4 . 0x10 2 °eV 
° o 


and 

a - 2xl0~ 2 8 cm 2 . 
o 

Then the source function for y-ray production may be written 


in the form 


q(E y ) = 4r 


2 a V< - < ) " r " 

E Y +m * /4 *y 


bb 


x dx dE I(E )a(E x)f(EjE ) 


(4) 


We assume all the pions to be produced at the average 
energy 




2 


m 2 + 2e'M 

TT p 

M' 2 + 2e 'M 
P P 


(5) 


Using equation (4) , equation (5) further reduces to 
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E 


E m 2 + 0 . 7M 2 

TT P 


; _E 

2 


1.7M 


“ V 5 


(6) 


so that the distribution function is approximated by 


f(EjE ) -6 < E n~ V 5) 


(7) 


Further specifying the differential cosmic-ray proton 

-r 

spectrum by a power-law of the form I(Ep) = K E equation (4) 

P P 


reduces to 


q (E ) — 8 TTn, , 0 K 
^ Y bb o p 


d V E J x ax[dE e'^-e /E /5) 

4 E Jo k /2) 


Y 


Max [(E q /5E ) ,2] 


( 8 ) 


8TTn,,c K (5E r ) / x (F 1) dx 
bb o p o / 

u 9 

0 


The solution to equation (8) may then be written in the 
simple form 


f Q 


q(E y ) = 1 


Y 


' M E i 


-r 


L Q y {E y /E 1 ) , E y *Ej 


( 9 ) 
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where E = E /10=4 . 8x10 9 eV and 
1 o 


d-r) 


Q = 8im. .a K r 
bb o p T 


E ^ = const. 


( 10 ) 


These photons travel, on the average, one mean free 
path A^E y ) before converting into an electron-positron pair 


1/A (E ) = (Tra 2 /3A) (kT/in ) 3 (m VkTE ) inf 0.117) (kTE v /m )] 
PP T e e y l i ^ j 


where A=ft/m c and a=l/B7=e 2 /«c (Gould and Schrader 1967) 
e 

The first-generation y-rays reaching the earth from 

distances up to A then have a flux given by 

PP 


I (E ) = q (E ) A (E )/4tt 
i ' y ^ y PP Y 


( 12 ) 


From equations (9) , (11) and (12) 


I. (E ) = L Lin(E v ) i 


f(E /E ) 
‘ Y i 


, E < E 
Y x l 


l Y 


Y 


4u 


{(E Y /E i ) _(r " 1) , E y >^ E i 


(13) 


where 


Q 

4tt 


8 X io' 25 k _ 3 

— = * cm s sr eV 


2 o r 

T(2.4X10 ) 


with 


(ID 


( 14 ) 
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H(E ) = 4.34/(lo gi 0 E y -15.97) 


( 15 ) 


and L = 1.0xl0 2 5 cm. 
i 

As discussed previously, photons undergoing pair con- 
version interactions after traveling distances vA^or greater, 
give up most of their energy to one electron. The fraction 
of energy given off to the other electron is given by 


where 




In 2 + 1/4 
*n(E /£)-! 


(16) 


£ = m Z /2e - 2.08 X 10 1 4 eV 
e o 


If the intergalactic magnetic field is weak enough 
so the Compton interactions with the microwave radiation 
dominate over synchrotron radiation as an energy loss mechanism, 
as discussed by Wdowczyk, et al. (1972) in the Klein- Mis hina 
high-energy limit (E^>>£) , the electron will transfer most 
of its energy to a microwave photon in a Compton interaction, 
converting it into a secondary ultrahigh energy photon. The 
difference in energy between the incoming electron and outgoing 
photon is 

E e - E y \ . 4/3 

E e / £n(E e /£)+l/2 



( 17 ) 
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In this manner cascade secondary photons can be 
generated producing a total flux greater than 1^ (E y ) . 
Wdowczyk, et al . (1971a, b; 1972) have made numerical calcu- 

lations of the total flux, I z (E y ) generated in the cascade. 
Here, we use a different approach to derive a simple approxi- 
mate solution for I z (E y ). The key to our approach lies 
in the fact that the total effect of the two-stage pair 
production-Compton process is to replace a photon of energy E 
by one of energy E y (l-fi ) (1- 6 C ) • Under the conditions 

S <<1 and <5 <<1, the two stage process can be treated as a 
pp c 

continuous energy-loss at a rate 


u(E y ) = dE y /dt 


(18) 


and I y can be approximated by a solution to the diffusion 

u 

equation 


9_ 

8E 


Y 


4i ui. 


- f Y q(V 


with the solution 


w = 


c f ( 

4ir |u(iy )| \J 


* A 

dE y q (Ey ) 


(19) 


(20) 


where, in equations (19) and (20), f y represents the mean 
portion of primary photons which are at any given time in the 
form of secondary photons rather than secondary electrons. 

This fraction is derived from the ratio of lifetime of photons 
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against pair-production x to that of electrons against 

Compton scattering x_ where, for E =E » £, 

^ 0 Y 

(x /x ) = (o /a )-> — t 1 } ^ — AZj. 

' pp 7 c ; E e =E y =E 1 c 7 pp' 2 £n(E/£)- 1 (21) 

Thus, 

' f Y " T pp/ <T PP +T c’ “ °' 37 f ° r lt>17 * ES1 ° 15 <22) 

From equations (16) , (17) and (21) , we find that the 
mean rate from the two-stage process is given approximately 
by 


u (E ) - 1.26 (E ) E c/A (E ) (23) 

y PP Y Y PP Y 


Substituting equations (11) and (16) in equation (23) and 
equations (22 and (23) in Equation (20) we find 


( 


W ” TiT 5< E yW ! 


r/(r-i) - (e^e^ ,e «e x 


(24) 


[l/(r-l)j (E^E,)"* 1 '"^^ 


where 

?(E y ) = ( log ! 0 E y -14 . 15 ) / (log j 0 E y -15.97) 


and L = 5. 9x10 2 5 cm. 

u 
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Equations (18) and (24) are valid in the energy region 
where the 2.7K radiation field is the dominant source of 
target photons and where the production of muon pairs is 
unimportant, i.e., for E y <10 19 eV. Wdowczyk et al. (1972) 
have determined that for 10 1 9 <E y <10 2 c eV X pp (E y )=10 25 cm* 
const. Thus, in this energy range, equations (13) and (24) are 
overestimates and have to be corrected for these additional 
absorption processes. Also, in equation (24) we have neglected 
the additional contribution of the lower energy electrons 
produced in the pair-production process. This correction 
becomes important at lower energies where, e.g. equation 
(24) underestimates I^E ) by a factor of ^2 at E y =E / 6 pp (E x } “ 

E /20=2.4xl0 18 eV. Estimates of the above corrections have 
been made to equations (24) and the resultant y-ray spectra 
obtained for the ultrahigh energy primary spectra given by 
Linsley, (1962), and Andrews, et al. (1971) are shown in 
figures (1) and (2). To the extent of their validity, the 
results of figures (1) and (2) confirm the numerical calcu- 
lations and conclusions first obtained by Wdowczyk, et al. 

(1972) for photon spectra in the energy region 10 1 8 <E y <10 2 °eV. 

The ef.fect of magnetic fields on the cascade process 
was also pointed out by Wdowczyk, et al. It becomes important 
when the electron energy loss rate due to synchrotron 
radiation in intergalactic space is greater than or 
equal to that due to Compton interactions, a condition fulfilled 


when 
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B. >> 1.3 X 10 -11 (E . /E ) gauss (25) 

ig 1 e 

A recent discussion of the possible intensity of 
intergalactic magnetic fields has been given by Brecher 
and Blumenthal (1970) . Using their arguments, but with 
the new value of the Hubble constant of 50 km/s/Mpc, it 
follows that if the general galactic magnetic field condensed 
out of the intergalactic field then B^ > 10 10 gauss. 

However, if the contrary is true and intergalactic fields 
result from expansion of galactic fields then B_^ a <10 8 gauss. 
A proper determination can only be made empirically and the 
empirical upper limit on the intergalactic field based on 
determination of the rotation measure for distant radio 
sources is 


B ■= £ 2 X lO -10 gauss (26) 

ig 

Thus, from equations (25) and (26) it follows that the 
cascade process may or may not be cut off at energies 
greater than or equal to 3x10 18 eV depending on the exact 
value of the average intergalactic magnetic field. 

However, the galactic magnetic field has a much more 
predictable effect on the ultrahigh energy electrons which 
make up 1-f =63% of the electromagnetic cascade. These 
electrons have to travel through at least half the thickness 

of the galactic disk *6.5x10 2 °cm spending *2x10 1 °s 
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in the galactic disk with a mean field of %3xl0 6 gauss. 

This is enough time to reduce the energy of any electron 
with E e >>10 15 eV to an energy of -10 15 eV (out of our region 
of interest) by synchrotron radiation losses. Thus, regardless 
of the strength of the intergalactic magnetic field, the 
galactic field strength ensures the elimination of all 
ultrahigh energy electrons from the cascade before reaching 
the earth. 

The electron and neutrino flux from the decay of 
charged pions produced in photomeson interactions will now 
be estimated. Again, the main contribution is dominated 
by the channel 

f n + tt + Prob.= 1/3 

y + p + A ( 1238 ) -> -! ^ (27) 

l p + tt° Prob.= 2/3 

and thus for rough estimates we take 

- 10 -28 cm 2 (28) 

The decay of the tt + meson results in a positron and 
three neutrinos all with an average energy v35MeV. The 
positrons represent an additional contribution to the 
electromagnetic cascade flux of the order of 10-20%. The 
neutrinos provide a source spectrum of the form similar to 
eequation (9) but with E i replaced by E^^E^/2 and with 



( 29 ) 
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so that at high energies 

Q<V 3 ( W~ r p 3 

^7 E =E =E * "5 <V 2 V‘ r “ 15 
The resultant neutrino flux is given by 


( 30 ) 


I (E ) = -r-^—R - 12 . 51 (E )f or E =E >2.4xl0 19 eV (31) 

v v 4 tt u P v p 

Where R =10 2 8 cm is the radius of the universe since the 
u 

universe is fully transparent to neutrinos. This possible 

large flux of high energy neutrinos has been suggested by 

Berezinsky and Zatsepin (196 9 a, b) to be a possible 

explanation for air-shower events of energy greater than 10 2 °eV 

if the neutrino interaction cross section becomes comparable 

to that of hadrons at these energies. Such an explanation 

is needed if these events are to be of extragalactic origin 

since the cosmic-ray proton spectrum is expected to cut off 

at E =10 2 °eV (Greisen. 1966,. Zatsepin and Kuzmin 1966,Stecker 1968). 
P 

Assuming then that the neutrino inelastic cross section 
rises linearly with energy and is given by (Krishnaswamy 
et al. , 1971; Bjorken and Paschos 1970 rPattison 1970) 


a 


- 10 


-38 


! GeV 3 3xlO' 2z a p>lnel E e . v (32) 


we find 


v 


P 


— at E- 2 .5x10 20 eV 
12.5 


( 33 ) 
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Since the extragalactic proton spectrum should drop 
off at 10 20 eV and the nutrino event rate va I , the event 

V V 

rates for neutrino induced and cosmic-ray induced showers 
may become comparable in the 10 2 °eV energy range but with 
a predicted flattening in the spectrum by one power at the 
higher energies. The transition energy may be somewhat 
lower depending on the exact values of r and and the 

role of cosmological neutrinos. (See discussion of Berezinsky 
and Zatsepin 1969b and Stecker 1971*). 


*In a previous paper (Stecker 1971) , the author considered 
the cosmological neutrino hypothesis unaware the Berezinsky 
and Zatsepin (1969b) had ruled out the same form of that 
hypothesis on the same grounds in a reexamination of their 
earlier paper (Berezinsky and Zatsepin 1969a. Those arguments 
do not apply to the non-cosmological neutrinos discussed 
here and by Berezinsky and Zatsepin. The critical test 
for the non-cosmological neutrino hypothesis lies in a 
flattening of the ultrahigh energy air-shower spectrum 
which may (Linsley 1962) or may not (Andrews, et al , 1971) 
exist. 


In conclusion, it appears that it may yet be quite some 
time before the uncertainties regarding the various possible 
components of the ultrahigh energy cosmic radiation are 
finally resolved, but the following points may be made: 
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1. The presence of photons making up a few per cent 
of the total flux at 10 1 9 eV would .imply a universal flux 
of both blackbody photons and ultrahigh energy cosmic-ray 
nucleons and also indicate that the intergalactic magnetic 
field has an average value less than <vlO 1 1 G. On the contrary, 
the absence of a significant 10 1 9 eV photon flux does not 

rule out the universal ultrahigh energy cosmic-ray hypothesis 
since it may only indicate the presence of intergalactic 
magnetic fields >10 1 1 G. 

2. Ultrahigh energy electrons are not expected to 
reach the earth because of synchrotron radiation losses in 
the galactic magnetic field. 

3. The establishment of evidence for a large ultrahigh 
energy neutrino flux would indicate the presence of both 
ultrahigh energy cosmic-rays and black-body radiation 
throughout the universe. The absence of neutrino events 
(with such events evidenced by a flattening in the air-shower 
energy spectrum) may only indicate that a does not rise 
linearly with up to the highest observed energies. 

4. The presence of a continuing straight primary 
spectrum above 10 2 0 eV would indicate a local origin for 
these particles either in the galaxy (Syrovatsky 1971, 

Stecker 1971) or in the local supercluster (Stecker 1968, 
Brecher and Burbidge 1971) , or (least likely) the absence of 
a universal blackbody radiation field. In the case of 
galactic origin, the ultrahigh energy cosmic-rays are expected 
to be heavy nuclei (Stecker 1971) . Some evidence for this 
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has been summarized previously (Stecker 1971) ; more recent 
evidence that the ultrahigh energy cosmic-rays may be heavy 
nuclei comes from a reexamination of the Volcano Ranch 
data (Linsley 1972, private communication) and a recent 
suggestion by Wdowczyk and Wolfendale (1972) based on the 
scaling hypothesis (Feynman 1969) . 



References 


Allcock, M.C. and Wdowczyk, J. 1972. II Nuovo Cimento 9B, 315 . 
Andrews, D., Edge, D. M. , Evans, A. C., Reid, R.J.O, Tennent, A 
1971, Proc. 12th Int'l. Conf. on Cosmic Rays . (Hobart) , 993 . 

Berezinsky, V.S. and Zatsepin, G.T. 1969a, Phys. Lett. 28B , 423 . 
Berezinsky, V.S. and Zatsepin, G.T. 1969b, Yad.Fiz. 10,1228. 

Translation: Sov. J. of Nuc . Phys. 10,696(1970) . 

Bjorken, J.D. and Paschos, E.A. 1970. Phys. Rev. Dl , 3151. 
Brecher, K. and Burbidge , G. R. 197 2 , Ap. J. 17 4 , 253. 

Brecher, K. , and Blumenthal, G. R. 1970, Astrophys . Lett. £,169 
Bonometto, S. A. 1971. II Nuovo Cimento Lett. 1, 677. 

Feynman, R. P. 1969. Phys. Rev. Lett. 23, 1415. 

Gould, R. J. , and Schreder, G. P. 1966. Phys. Rev. Lett. 


16, 252. 


, 1967, Phys. Rev. 155, 1404. 

Jelly, J. V. 1966. Phys. Rev. Lett 1£, 479. 
Linsley, J. 1963, Phys. Rev. Lett. £, 126. 
Krishnaswamy , M. R. , Menon, M. G. K. Narasimham, 
Hinotani, K. , Ito, N., Miyake, S., Osborne, 
Parsons, A. J. , and Wolfendale, A. W. 1971. 
Soc. Lond. A3 2 3 , 489. 

Pattison, J. B. M. 1970 Acta Phys. Sci. Hung. 29 
Stecker, F. W. 1968. Phys. Rev. Lett. 1£, 1016. 

1969. Phys. Rev. 180, 1264. 
1971. Nature Phys. Sci. 234 ,28 


V. S. , 

J • L . , 

Proc . Roy . 

Suppl . 4, 337. 


Syrovatsky, S. I. 1971 Comments on Astrophys. and Space 



Wdowczyk, J. , Tkaczyk, W., Adcock, C. and Wolfendale, A. W 


1971a. J. Phys. A4, L37 . 

1971b. Proc. 12th Int'l. Conf. on Cosmic Rays 

Hobart) 965. 

Wdowczyk, J. , Tkaczyk, W. , and Wolfendale, A. W. 1972 
J. Phys. A, 5_,1419. 

Wdowczyk, J. and Wolfendale, A. W. 1972. Nature Phys. Sci 
236, 29. 




Figure Captions 


Figure 1 


I 


Figure 2 


Cosmic-Ray Photon Cascade spectrum (I ) and first 

generation spectrum (Ii) for a primary spectrum 

(I shown) of 
P 

p (E p ) = 6 . 4xl0~ 4 °(E p /10 2 0 ) ~ 2 * 6 (cm* s • sr • eV) ” 1 
(Linsley 1962) . 

Cosmic-ray Photon Cascade spectrum (I ) and first 

generation spectrum ( I i ) for a primary spectrum 

(I shown) of 
P 

Cp(Ep) = 2.4xl0 - *°fE p /10 2 0 )" 3 - 24 (cm's.sr.eV)' 1 


(Andrews, et al . 1971) . 
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